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Abstract 

Using a generalisation of an Onsager type approach, we are able to predict a dielectric permittiv- 
ity profile of an inhomogeneous dipolar fluid in the presence of a dielectric interface. The reaction 
and cavity fields are calculated semi-analytically using bispherical coordinates. An asymptotic 
expression for the local permittivity is derived as a function of distance from the interface. 
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I. INTRODUCTION 



Since Onsager ^ and Kirkwood [3| it is well known that the macroscopic dielectric 
constant (or permittivity) e of a material of polar, polarizable molecules is related to the 
fluctuations of the total dipole moment of a spatially homogeneous sample [^. While the 



relation is well understood for infinite or periodic Q bulk materials, much less is known 
about the static dielectric properties of confined polar fluids, or polar fluids near interfaces. 
The generally accepted view is that near interfaces, the dipole moment fluctuations are 
restricted, e.g. when dipole moments of individual molecules tend to orient themselves pref- 
erentially parallel or perpendicular to an interface, and that the resulting local permittivity 
is reduced relative to its bulk value. This is particularly important when considering water 
molecules near electrodes, biological macromolecules (e.g. proteins or DNA) or membranes 
in electrochemical and biophysical applications. On purely phenomenological grounds it is 
often assumed that the dielectric response of a polar solvent in such situations may be de- 
scribed by a distance-dependent permittivity, which is then incorporated in expressions for 
the Coulomb interaction energy between ions in solution 0. 

The present paper attempts a first step beyond such purely ad hoc procedures, in terms 
of a simplified molecular picture of a polar fluid near an interface, by generalising Onsager's 
well-known cavity model for bulk dielectric properties to an inhomogeneous situation. The 
objective is to determine a permittivity "profile" which varies with the distance from a 
planar interface separating a continuous medium of given permittivity ei from the medium 
made up of dipolar molecules. 

The paper is organised as follows: The generalised Onsager cavity model is described in 
section |n] A self-consistent set of equations determining the permittivity profile is derived 
in section Hm In section HVl the asymptotic behaviour of the dielectric permittivity far away 
from the interface is examined. The solution to the full electrostatic problem leading to 
a locally self-consistent permittivity profile is sketched in section |^ The results will be 
discussed in the final section. 
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II. ONSAGER MODEL NEAR A PLANAR INTERFACE 



Consider two dielectric media separated by a planar interface, as shown in figure The 
horizontal interface is placed at 2; = 0. Below it {z < 0) extends a dielectric medium 
considered to be a homogeneous continuum of constant permittivity ei. Above the interface 
{z > 0) extends a polar fluid made up of identical spherical molecules carrying a dipole 
moment /x of fixed magnitude fi = The interface breaks the rotational invariance of the 
fiuid, so that its dielectric response is a priori described by a permittivity tensor €2 which, 
by symmetry, is diagonal and of the form 

0^ 



£2 



el 



(1) 

yO 

where e| denotes the identical xx and yy components parallel to the interface, while 
denotes the vertical zz component. 

In view of an approximate calculation of e| and e^, one assumes, following Onsager Q|, 
that a test molecule is placed at the centre of a spherical cavity of radius R = (3f/47r)^/^, 
where v = V/N is the volume per molecule. Since we are interested in permittivity profiles, 
the centre of the cavity is placed at a vertical distance z from the planar interface, and the 
dielectric response will be determined as a function of z. The cavity is surrounded by a 
dielectric continuum characterised by the dielectric tensor £2, of yet unknown components 
e| and e^. Following Onsager's mean field approach 1], the dipole /x inside the cavity is 
subjected to a cavity field E^, induced by some external field taking the value E'^^* far above 
the planar interface (i.e. at z +00). In addition, the dipole polarises the surrounding 
dielectric media I and II, giving rise to polarisation charges on the surface of the cavity and 
on the planar interface. These polarisation charges generate a reaction field Er acting on the 
dipole inside the cavity. Both fields E^; and E^ are calculated from the laws of macroscopic 
electrostatics. The mean dipole moment inside the cavity is then calculated by averaging 
over orientations with the appropriate Boltzmann weight determined by the coupling of 
the dipole to the total electric field E = Ec + E^ inside the cavity. The dielectric tensor 
components e| and are finally determined by identifying the calculated mean dipole with 
the predictions of macroscopic electrostatics for parallel and orthogonal orientations of the 
external field. 
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I 

1 

FIG. 1: Geometry of the Onsager model. The test dipole is in a spherical cavity with centre 
a distance z above the dielectric interface, below which extends a continuum of permittivity ei. 
The fluid in region II has a local permittivity tensor £2. Far above the wall the electric field is the 
external field E''^*. The polar angle of dipole with respect to the z axis is denoted by 

The calculation sketched above is quite straightforward in the case of isotropic, bulk 
dielectrics, but, as will become clear below, it is technically much more involved in the 
presence of an interface between two different dielectric media. 

In the isotropic bulk, where e| = = e2, the Onsager mean field theory leads to the 
classic result Q 

(62-1) (262 + 1) ^ _p^^ 

valid for non-polarizable, polar molecules, p = 1/v being the number density of these 
molecules and ksT the thermal energy; 60 is the vacuum permittivity. The predictions 
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of (j2I), which does not account for correlations between the orientations of neighbouring 
dipoles (as embodied in the Kirkwood (^i^-factor, see eq. ()A2|) ). considerably underestimate 
the value of €2 for dense, highly polar liquids, and the same is hence to be expected for the 
results near the planar interface. 



III. SELF-CONSISTENT EQUATION FOR THE PERMITTIVITY PROFILE 

In the presence of the planar interface (cf. Fig. Q), the polarisation vector P in region II 
depends on the vertical coordinate z. If £(2;) is the local electric field, the local susceptibility 
and dielectric tensors X2 ^-nd can be defined as usual by the linear relation 

V{z) = eoX2i^) ■ E(^) = eo(e2(^) - I) ■ E(z), (3) 

assuming a purely local response. 

As explained in the previous section, the local field acting on the dipole ^t inside the 
cavity is the sum of the cavity and reaction fields. Due to the linear nature of electrostatics, 
the reaction field depends linearly on the dipole moment 

= Mr ■ fi. (4) 

The tensor Mr depends on ei and £2 and on the geometric parameters R and z. Due to 
the cylindrical symmetry of the problem, it is clear that Mr must be diagonal. A vertical 
dipole (parallel to the symmetry axis Oz) can only induce a vertical reaction field. For an 
arbitrary orientation of /x, the reaction field must change sign under the inversion of the 
dipole i-i — > —/X. In the case of a horizontal dipole, this inversion is equivalent to a rotation 
by vr around the vertical symmetry axis, implying that the reaction field has no vertical 
component. The horizontal component must moreover be parallel to the dipole, since the 
system is invariant under refiections by the plane containing /x and the vertical symmetry 
axis. The tensor Mr is hence of the same diagonal form as the permittivity tensor, with 
one vertical component = Mrzz-, and two horizontal components m\ = Mjixx = ^Ryy 
Similarly, the cavity field is proportional to the external field 

Ec = Mc ■ E^^\ (5) 

where the tensor Mc can only depend on ei, £2, and on the geometric parameters R and z. 
The same symmetry arguments as adove show that the tensor M^ is diagonal. 
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If the dipole interact with the cavity field via the energy — /x ■ Ec, the interaction energy 
with the reaction field is given by — /x-E/j/2. Indeed, increasing gradually the dipole moment 
from zero to its final value, we obtain a factor 1/2 when we integrate the infinitesimal 
interaction energy increment — E/j ■ because of the linear relation (^. 

The induced dipole moment per unit volume (or polarisation) is, within the cavity model. 



where the integrations are over all orientations of the test dipole. For small external field, 
the polarisation will be proportional to the field, and we can expand the exponential to first 
order in the reaction field. Taking into account the diagonal nature of the cavity field tensor 
Mc, we obtain in linear response 

J tiexp{f3ER ■ ^i/2}{l + f3Ec ■ ti)dQ 
J exp{/3En ■ ^/2}(1 + (3Ec ■ ti)dQ 
J fxexpipEn- fx/2}l3Ec- fidn 
J exp{(3En ■ ^i/2}dQ 
Jfifiexp{(3^ ■ Mr ■ fi/2}dn 
Jexp{(3fl■Mn■^l/2}dn ^ 

eoXr-E^^'- (7) 



In going from the first to the second line of ((Tj), integrals involving an odd function of fi 
vanish. In the last expression we have defined the dielectric susceptibihty with respect to 
the external field 

ext ^ /3p/i^ / AAexpj/j/x ■ Mr ■ ^il2}dn 

eo jexp{(5^i-MR- ^l/2}dVL ' ^' ^' 



The components of this tensor are 

xT = ^-{1-^r)mI xr = v^ijM^, (9) 



with y = Ppfi'^/eo and 



1 e""" 



a=^/^i^{Ml-M^)/2. (10) 



20^ A/vraErf(a) 

(Erf denotes the error function). The external susceptibility (jHl) must now be related to the 
local susceptibility X2 defined in (jS)). 



In the present case of a planar interface, it is easy to find the macroscopic relation between 
the polarisation P{z) and the external field. For a given field far away from the wall and 
a given permittivity profile, the electric field at every point is known from the macroscopic 
equations. When the external field is horizontal, it will stay constant throughout space, 
since the permittivity is only varying perpendicularly to the field. Therefore the field felt 
by the test dipole in the macroscopic picture is simply E*^""*, leading to the polarisation 
P = eo(e|(2;) — 1)E^^*. In the case of a vertical field, however, the displacement field 

D(^) = eoE(z) + P{z) = eo 62 • E(^) (11) 

will remain constant. Therefore, at distance z from the interface, the electric field is 
E''^*e^(oo)/e^(;z), leading to the polarisation P = eo(e^(oo) - ^^)E'=^*. The permittivity 
far away from the wall is simply the bulk permittivity. 

Comparing the above expression of the polarisation with the result (jH]) obtained from 
Onsager's model gives 

X^2l = 4W-l = y(l-Aii)Mll (12a) 

X?1 = 4(oo)-^ = |a^M^. (12b) 

ef{z) 2 

The tensors Mc and M/j [which determine A/j via (fTUI) ] are to be obtained by solving 
the associated electrostatic problems. Since the permittivity tensor £2 (-2) itself enters these 
problems, equations (jl2ap and ()12b|) are self-consistent equations for this permittivity profile. 

Before determining the cavity and reaction field tensors in section V, we obtain first the 
asymptotic behaviour of the permittivity. 



IV. ASYMPTOTIC BEHAVIOUR OF THE DIELECTRIC PERMITTIVITY 

It is instructive to obtain an approximate estimate of the corrections to the reaction 
and cavity fields, valid for a test dipole far away from the interface, such that z ^ R. 
The dielectric tensor will then practically reduce to the bulk value e2l. The polarisation 
of the cavity surface by the test dipole in the absence of an external field is such that 
the latter acts outside the cavity as an effective dipole moment of strength The 
cavity model only yields an approximate estimate of the effective dipole moment. An exact 
calculation, sketched in Appendix |X1 yields a general expression for the screened (effective) 
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dipole moment which reduces to the prediction of the cavity model if the Onsager value (j21) 
for €2 is substituted in equation ()A6|) . 

The additional electric field generated by the polarisation of the planar interface aX z = 
by this effective dipole may be described as being due to an image dipole of strength 
situated at —z below the interface. Since z/R^ 1, the resulting field may be 
regarded as homogeneous over a distance equal to the diameter 2R of the cavity. The reaction 
field inside the cavity will be that outside the cavity, amplified by a factor 3e2/(2e2 + 1). 
The field depends of course on the orientation of the image dipole which itself is determined 
by the orientation of the test dipole. If the test dipole has both vertical and horizontal 
components, /x-*- and /x" , the reaction field is therefore 

2 



E 



1 



R 



47rer 



^62 — 1 ^ 

'262 + 1^" V2e2 + 1 



''2 ~ ^1 1 r II , o J-i 
£2 — : — — [m" +2/^ J 



(13) 



■e2 + ei 82;3 

and hence is no longer coUinear with /j,. The electrostatic interaction energy of the dipole 
with the reaction field — ^t ■ E/j therefore depends on the orientation of the molecule. In 
figure 121 it can be seen that the orientation dependence of the energy is negligible even in 
the immediate vicinity of the wall. 

In the presence of an external field two parts are contributing to the effective dipole 
moment. The test molecule aligns on average with the cavity field, leading to a permanent 
dipole, which appears screened outside the cavity. A second contribution originates from 
the polarisation of the cavity surface by the external field. Both terms lead to an electric 
field outside the cavity as due to an effective dipole moment 



.cff 



11^ 



3e2 



2e2 + 13A;BT2e2 + l 



47ref 



62-1 
'2eo + l 



E' 



ext 



(14) 



The corresponding image dipole /i™^ 



62-61 ,.efF 
e2+ei 



/Li™ induces an approximately homogeneous field 



E™^ near the cavity, which is the first order correction of the local electric field E "'^ = E™^ + 
E*^^* over the external field. From the local field the local polarization can be determined to 
first order in z~^ 



9eol/ 



2e2 + l" 



£2 



262 + 1 



E 



cxt 



J^ext _|_ 2^6xt 



± €2 - ei 



327760^^ 62 + ei 



3e2/3/i' 
L(2e2 + i: 



47reo-R' 



3 ^2 



262 + 1 



(15) 
(16) 
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FIG. 2: Potential energy —fi ■ divided by the thermal energy ksT as a function of the 

azimuthal angle "i? due to the interaction of the dipole with the reaction field as predicted by the 
Onsager model (solid line) and as given by the asymptotic expansion ((T^ (dashed line) . The dipole 
is located in a distance 2R above the wall. Dipole moment /i = 1.86D = 6.2 • 10^^'^Cm and number 
density p = 3.346 • lO^'^m"^ are the experimental values measured for water at T = 293K. This 
leads in the bulk to an Onsager value of 62 = 18.4. The material below the interface is vacuum 
ei = 1. The constant electrostatic interaction energy which is also present in the bulk has been 
subtracted. The amplitude of the orientation dependent energy is less than COlfc^T even in the 
immediate vicinity of the wall. 



After extracting xT can find the dielectric permittivities to first order in z as 



£2(00) 



€2-1 €2 - ei 



327reoz3 62 + ei [(262 + 1)^ 



3e2/5/i2 



3 ^2 



269 



- £2(00) 



62(62 - 1) £2 - ei 
IGneoZ^ 62 + Ci 



(2e2 + 1)^ 



47r6o-R' 



; ^2 - 1 

262 + 1 



(17) 
(18) 
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All quantities appearing on the r.h.s in these expressions are bulk values. 

V. SOLUTION OF THE ELECTROSTATIC PROBLEM 

To find the tensors M(7 and we have to solve two purely electrostatic problems. One 
involves finding the reaction field induced by the polarised material at the centre of the 
cavity. In this case we consider a dipole at the centre of the spherical cavity, surrounded 
by dielectric material of permittivity 62 at a distance z above a medium with isotropic 
permittivity ei. The electric field must decay to zero as 2; — 00. To find the cavity field on 
the other hand, we consider an empty sphere, and an external field which goes to E''^* far 
from the cavity. In both problems the electrostatic equations are the same; they differ only 
in the boundary conditions: for the cavity field we require the electric field to be finite within 
the sphere and approach a prescribed value at infinity; for the reaction field we require a 
dipolar singularity at the centre of the sphere and the electric field vanishes far away from 
the interface. 

In the geometry considered here there are two dielectric interfaces: One is a planar 
interface where the dielectric permittivity jumps from t\ to e^. The other is a sphere where 
the permittivity jumps from to 1. These interfaces divide space into three regions as 
shown in figure ^ Within each region the electrostatic potential satisfies the relation 



With a spatially varying tensor permittivity €2(2;), eq. (jl9|) cannot be solved analytically. 
Instead of computing a numerical solution, we shall determine the cavity and reaction fields 
for the simpler case of a constant and scalar permittivity €2. These results will then be 
used to determine a permittivity profile from the Onsager equations ()12a|) and ()12b|) . These 
equations will be solved in a locally self-consistent manner: when the cavity is at a distance 
z from the interface, we shall use the cavity and reaction field tensors corresponding to the 
case of a homogeneous and isotropic medium II of dielectric constant = Tr(e2(z))/3 = 
[2e|(z) + e^(2;)]/3. This scheme provides a reasonnable approximation to e{z), since the 
deviations of e{z) from the bulk behaviour €2! prove to be small for z > 2R (see section VI). 
In each regions, equation (fT^ now reduces to Laplace's equation 



V ■ (e ■ V^j) = 0. 



(19) 



Aip = 0. 



(20) 
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At the dielectric interfaces we impose the usual boundary conditions: The tangential com- 
ponent of the electric field E = —'Vip must be identical on both sides of the interface: 



(21) 



t,i 



tji 



This condition is automatically satisfied if the potential is continuous across the interface: 



(22) 



Likewise the normal part of the displacement field D = eoe2E must be identical on both 
sides of the boundary. 



62 Vil) 



nj 



(23) 



njl 



The electrostatic problem is best solved in bispherical coordinates. Some properties of 
this system are summarised in Appendix ^ This coordinate system is well suited for our 
problem because the two interfaces turn out to be surfaces of constant coordinate rj. A set of 
these surfaces is sketched in figure 0] Furthermore Laplace's equation is separable in these 
coordinates. We therefore proceed to express the yet unknown electrostatic potential in each 
region as a series of fundamental solutions of Laplace's equation with unknown expansion 
coefficients. The latter will be determined by the boundary conditions. In bispherical coor- 
dinates the component of the displacement field normal to the surface of constant coordinate 
T] is given by |lO| 

coshr] 



Dn 



COS i!) dip 
or] 



(24) 



We fix the parameter a in the bispherical coordinates to a = Ry^ z'^/R'^ — 1 and define 
rjo = \n{z/R+ {z/ RY — 1). With these choices the planar interface corresponds to the 
surface = 0, and the spherical cavity to the surface rj = rjo. The centre of the sphere has 
coordinates t]i = 2r]Q, = 0, and 0i = 0. The solutions of Laplace's equation in bispherical 
coordinates are of the form 



(25) 



n=0 m=—n 

with the fundamental solutions 



^iniv.^.'P) = Vcoshr/ - cos7?e±("+i/2)r,y™^^^ 



(26) 
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/\ z 

r| = 0.5 




r| =0 

.^Ti=-0.5 

FIG. 3: Surfaces of constant rj in bispherical coordinates. A cut through the cartesian xz-plane 
is shown. In three dimensions the surfaces are spheres. As 77 ^ 00, the centres of the spheres 
approach the point z = a. 

Here the are spherical harmonics and the A^n and Bmn are constant expansion coef- 
ficients. The ipmn diverge at r] ±00, i.e. x = y = 0, z = ±a. Because the potential 
in region I is regular, only the ip^,^ can contribute to the expansion. Likewise the expan- 
sions of the cavity and reaction field inside the cavity only contain terms. However, 
the field produced by the dipole itself contains contributions from the divergent ip:^,^ (see 
Appendix • This is due to the singularity of the potential at the centre of the sphere. We 
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therefore expand the potential in each of the regions in the following series: 

oo n 

^I = Y.Y. ^rnn^^n (27) 
71=0 m=—n 
oo n 

i^n = J2Yl [Bmni^^n + Cmni^i^] (28) 
n=0 m=—n 
oo n 

i'lll = XI 5Z [^rnn^^n + ^mn^^^n] > (29) 
n=0 m=—n 

bearing in mind that the Emn are the known expansion coefficients of the dipolar potential 
given by equations ()B9|) and ()B10|) . When calculating the reaction field, the Emn vanish. 
In this case the expansion coefficients of the external field, given by equations ()B11|) and 



()B12|1 . have to be included in the Bmn- The expansion coefficients are determined by the 
boundary conditions. Condition ()22j) . i.e. continuity of the electric potential, directly leads 

to 



oo n 



m. 

n=0 m=—n 

oo n 



oo n 



»?=0 71=0 m=—n 



77=0 



n=0 m=—n 



(30) 
(31) 



v=vo 



V=VO 71=0 m=-n 

for all and 0. Multiplying these equations with Y^^*{{}, 0) / v^co3ir7--"cos^ and integrating 
over sin-t^d-^d^), one is able to separate the terms in the sums, with the result: 



Amn — Bmn + Cmn 



(32) 
(33) 



The second condition, i.e. continuity of the normal component of the displacement field, 
eq. (|2S1), leads to the equations 



oo n 



^1 / / ^ ^mn 
71=0 m,=—n 
oo n 



£2 



EE 

71=0 m=—n 



drj 



ri=0 



oo 71. 

^^EE 

n=0 m=—n 

dip; 



" y mn i rmn 

-Dmn ^^^^ ~r ^5 



9?7 



drj 



in o 

or] 

oo n 

EE 



r?=»70 71=0 m=-n 



Dr. 



drj 

di/j- 



dr] 



77=0 



+ E„ 



Or] 



Projecting once more the sums, we find 



^iAmn — ^2{Cmn — Bmn)- 



(34) 



v=vo 
(35) 



(36) 
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Equations (jH^ . and can be used to express Amn, Bmn, and Cmn in term of the 
unknown -Dmn- 



A„ 

Br, 



2e2 



£2 + ei 
£2 - ei 



^ (D + E e^^"+^^*) 



with 



£2 + ei 

C = ^ (D + E e(^"+^)''«' 



1_ _ £2-61 ^ ^(2«+l),70_ 



€2 + ei 

Equation now reduces to 



(37) 
(38) 
(39) 

(40) 



n=0 m=— n 



+ n cosh ?7o 



(n+l/2)»?o 



0% 



n cosh TjQ 



— [(e2i?mn — Dmn)e 



-(n+l/2)»?o 



+ (E^„ - e2C^n)e("+^/2)*] (n + 1/2) cos^^jFr 0) = 0. (41) 

We multiply this equation with Y^* {{}, 0) and integrate over sin 'dd-ddcf). Using the product 
rule for spherical harmonics (see eq. (A. 26) of ^]), we find 



-rio 



+ N cosh ?7o 



[e2BMN - Dmn] e-^^+i/^)"" 
+ [Emn - e^CMN] e(^+^/2)''° 

OD n 

J2 E [(e2S™n-n„„)e-("+i/2)''° + (i?^„-e2C^.)e(""''/'^*] 



— + cosh ?7o 



n=0 'm=—n 



2n + l 



X {n + 1/2) ^^ ^^ C{nlN; 000)C{nlN; mOM), (42) 

where C{nlN]mOM) are Clebsch-Gordan coefficients. All terms with m ^ M vanish, as 
well as all terms with ri 7^ A^, A^ ± 1. The non-zero coefficients are 



C(nl(n + 1); mOm) 
C(nln; mOm) 



{n — m + l){n + m + 1] 



(2n + l)(n + l) 



m 



^yn{n + 1) 



C{nl{n — 1); mOm) = — ^ 



' {n — m){n + m) 
n{2n + 1) 



(43) 
(44) 

(45) 
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Equation therefore reduces to a three term recursion couphng the Dmn, Dm{n-i), and 
Dm{n+i)- When only a finite number of terms are taken into account, the resuhing tridiagonal 
hnear equation can be solved with standard methods. From the D^n the electric field inside 
the cavity can be obtained using equation ()B5|) . The series Yl'^=oJ2m=-n^rnn4'mniv^^^(p) 
converges very quickly at the centre of the sphere, so that a truncation after 20 terms leads 
to an accurate value for the field. 

After having calculated the reaction and cavity fields using a given dielectric permittiv- 
ity €2, equations ()12a|) and ()12b|) are employed to give new estimates for the permittivity 
tensor. A locally self-consistent permittivity can thus be found by an iterative procedure: 
We start by setting 62 to the isotropic bulk permittivity as determined by the bulk Onsager 
relation (j21). We then solve the electrostatic equations as outlined above, leading to xT^- 
Equations ()12a|) . ()12b|) are then used to obtain a new prediction of €2 (-2). The new isotropic 
permittivity €2 = Tr(e2(2;))/3 is then used as a starting value of the next iteration. In prac- 
tice this scheme converges very quickly. The resulting permittivity profile is shown in figure 

m 

VI. DISCUSSION 

By generalising Onsager's model for the dielectric permittivity in bulk fiuids we were able 
to predict a permittivity profile of a dipolar fiuid near a wall. A number of approximations 
have been made; some of which are of fundamental nature, and some are made to render 
the calculations more tractable. A fundamental weakness of the Onsager approach is the 
complete neglect of orientational correlations between neighbouring molecules. This leads 
in the bulk to a significant underestimation of the permittivity. In the vicinity of a dielectric 
interface the fiuid becomes birefringent. In our calculations we have neglected the anisotropy 
of the permittivity, to simplify the calculations. It is possible to take the full tensorial nature 
of the permittivity into account, which would lead to a five term recurrence formula instead 
of the three term recursion in equation However, the results (figure Hj) show that the 

fluid stays nearly isotropic except in the immediate vicinity (i.e. a few molecular distances) 
of the interface. Therefore the increased effort required to treat the anisotropy exactly 
does not seem to be worthwhile. In order to derive a semi-analytical expression for the 
permittivity, an inconsistency has to be allowed for in treating the model. For each position 
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FIG. 4: Dielectric permittivity profile as calculated self consistently. Dipole moment and density 
of the fluid are the values measured for water (see caption Fig. ^ . The material below the interface 
is vacuum. The inset shows that the dielectric permittivities approach the values predicted by the 
asymptotic expansion ((T7|) very slowly. 

of the cavity the permittivity of the whole fluid is assumed to be constant. The calculated 
permittivities, on the other hand, depend on the distance of the cavity from the interface, 
leading to a permittivity profile. This inconsistency can be overcome by solving the three 
dimensional Laplace equation numerically starting from a given permittivity profile. Again, 
the results of our model show bulk behaviour nearly everywhere except in the immediate 
vicinity of the interface, so that the complicated fully consistent solution does not seem to 
be worthwhile. Far away from the wall we were able to derive a first order correction to 
the bulk behaviour. We are presently investigating how to incorporate the orientational 
correlations between neighbouring dipoles into the theory. 
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APPENDIX A: EXACT EXPRESSION FOR THE EFFECTIVE DIPOLE MO- 
MENT 

Within a microscopic description of the dielectric medium modelled as a classical dipolar 
fluid ^j, we can obtain an exact expression for the effective dipole moment /-i'^^ of a given 
molecule inside the fluid. We assume the non polarizable molecules to have a permanent 
dipole moment /x, and we denote the position and the orientation of the ith molecule by 
i = (rj, Hj), i = 1, ...,N. The molecules interact via of the dipolar pair potential 

t;dip(l,2) = (/xi- Vi)(/X2 ■ Vs)^ r, |ri-r2|>0, (Al) 

|ri — 

and a short-ranged interaction fsr(l, 2). The short-ranged interaction (e.g. a Lennard- Jones 
potential) is quite arbitrary, except for the fact that it must decay faster than |ri — r2|^^ 
at large separations (it may thus include higher order multipole interactions between the 
molecules) . 

The dielectric constant e of this fluid, at inverse temperature (3 and number density 
p = N/V, is given by the Kirkwood formula 3| 

+ =yg^, y = 4nf3pf,y9, (A2) 

where the factor g^^ = 1 + p J d^ri f dQi f 6.^2 (Ai ■ ^2)^(15 2)/(47r)^ is a measure of short- 
ranged angular correlations between the molecules. The Onsager approximation to e, equa- 
tion (j21), corresponds to neglecting these correlations, i.e. setting g^^ = 1. We recall moreover 
that the Ursell function h{l, 2) of the infinite system decays asymptotically like a screened 
dipolar potential 



h{l, 2) ~ i--^(-/3t;dip(l, 2)), In - ral ^ 00. (A3) 



To evaluate n^^, we introduce p(l|2) = (p/47r)(l + h{l,2)), the density of molecules at ri 
with orientation fi^ when there is a dipole /Xg at r2. By definition, the average polarization 
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of the fluid around the flxed molecular dipole ^2 is 

P(ri) = J dn^ Ml p(l|2) = £-Jdni Ml h{l, 2). (A4) 

This polarization cloud carries a total dipole moment = J r (— V ■ P(r))d'^r. Ac- 

cording to ()A4|) and ()A3|1 . P(r) decays like |r|~^ at large distances and is therefore at the 
borderline of integrability. The integral deflning m'^'""*^ is nevertheless convergent, thanks 
to the harmonicity of the Coulomb potential: —V ■ P(r) decays rapidly since A(l/|r|) = 
when |r| > 0. An integration by parts gives 



Cloudy lim 



I P(r)d3r- I r(P(r) -dS) 

Jv JdV 



(A5) 



for any volume V. For a spherical volume of radius -R 00, we flnd using flA4j) and ()A2jl 
that the volume integral of the polarization is ((?k ~" 1)M2- The surface integral can be 
evaluated using ()A3|1 . with the result /X2 2(e — 1)^/ {9ye). This shows that the effective dipole 
moment is given by the simple result 

= M + A*'^""'^ = (A6) 
3ye 

APPENDIX B: PROPERTIES OF THE BISPHERICAL COORDINATES 



The bispherical coordinates [r], -t}, 0) are given by 



a sin -(9 cos . 
^ = — u q (^1) 

cosh 7] — COS V 



a sin sin 
cosh T] — cos {} 
a sinh t] 



(B2) 
(B3) 



cosh r] — cos -i? 

The surfaces of constant f] = rjo are spheres of radius R = a/\ sinhr^ol with centre at the 
Cartesian coordinates x = 0,y = 0, z = acothr^o- The bispherical coordinates of the centre 
(where the dipole is located) are r]i = 2r]Q, {}i = 0, and 0i = 0. 

Given the coordinate system ()B1|) one can work out the unit vectors in the 7], and 
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-direction as 



1 



cosh T] — COS ^ 



1 

cosh 1] — cos 'd 



^sin sinh r] cos (f^ 

sin sinh 77 sin q 
yl — cosh 7] cos ^ I 

^(cosh 1] cos -i? — 1) cos 



(B4) 



The gradient is 

, cosh ri — cos d 

= 



(cosh ?7 cos — 1) sin i 
sinh sin ■(9 



d'\\) dil) 1 9-?/' 

ar] ov sm -i/ 



^ - sin (h^ 



COS( 





(B5) 



The Laplace equation V^^ = is separable in the coordinate system and its solutions are of 
the form given in eq. ()25|) Moreover, the Green's function, i.e. the solution of the equation 
V2G(r) = -5(r-ri), is 

G(?7,'(9,0) = - A/cosh?] — cosT^-v/ coshr^i — cos-i^i- 
a 



EE 

n=0 m=— n 
ly y 

a ^ ^ 2n + 1 

n=0 m=—n 



2n + l 
Air 



(B6) 
(B7) 



(B8) 



The last line of the expansion is only valid for rji > rj. Thus we can derive the electrostatic 
potential of an elementary charge located at any point. From there we obtain the potential 
of a dipole located at the centre of the sphere r] = rjo. 

The potential for a dipole pointing in 2;-direction is obtained by taking the derivative of 
the Green's function 



V'dip('7)'*^)0) _ dG _ coshr/i - cos?9i dG 
dzi a dr] 



y/coshrji — 1 



E 

n=0 



An 
2n + l 



,-(n+l/2)r,i 



sinh rji 



- {n + 1/2) (cosh 771 - 1) 



(B9) 
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Similarly we find for a horizontal dipole in the x-direction 



r\-Jv,'^%<P) _ dG _ cosh Tji — COS 1?! dG 
H dxi a d-d 

n=l 



(BIO) 



Prom the generating function of the Legendre Polynomials we can find the expansion of the 
vertical external electric field with amplitude E^^^ valid in the region > 

oo 

i^Uv, ^, 0) = -Vs^aEr E \/2^^on(77, i^, 0). (BU) 

n=0 

Por the horizontal external field in x-direction we obtain 



4.t{v, ^, 4>) = 2v^aEf * V ^^^^ (^rn - ^-^n) ■ (B12) 

n=0 
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For a tensor permittivity 62, the factor 62 in this equation would be replaced by e|, and there 
would be an additional term -^^^^^^fS" V- - '/ 9^ 



(1 — cosh ?7 cos 1?) ^ + sinh?7sin^|^ 
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